The Magnetism of Li doped La 2 Cu04: the antiferromagnetic spin-shard state 
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We study the dynamics of a single hole in Li and Sr doped L^CuCU and its extension to a 
finite hole concentration. We compare the physics of La2- a; Sr :E Cu04 and L^Cui-^Li^CU and 
explain why these systems are remarkably different. We demonstrate that holes in I^Cui-^Li^C^ 
are always localized and that there is a critical concentration, x c « 0.03, above which the holes 
break the global antiferromagnetic state into an array of weakly coupled antiferromagnetic clusters 
(antiferromagnetic shards). We show that the spin-shard state provides a description of the magnetic 
and electric properties of La2Cui_^Li I 04. Two experiments that can test our theory are proposed. 

PACS numbers: 74.72.Dn, 75.10.Jm, 75.30.Fv, 75.50.Ee 
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Introduction. The sharp contrast in the physical prop- 
erties of Sr and Zn doping in La2CuC>4 is not surpris- 
ing since Sr brings mobile holes into the C11O2 planes 
while Zn does not. At the same time the remarkable 
differences between Li and Sr doping of La2Cu04 are 
somewhat surprising since both Sr and Li ions bring mo- 
bile holes into CuC>2 planes. It is by now well estab- 
lished that La2_^Sr a; Cu04 (LSCO) and La2Cui_ T Li 2; 04 
(LCLiO) have very different physical properties: LSCO 
has an insulator-superconductor transition at x « 0.055 
1[ while LCLiO remains an insulator at all dopings 
2j. Elastic and inelastic neutron scattering in LSCO 
at x > 0.02 reveals incommensurate magnetic peaks 
El 01 El while neutron scattering in LCLiO is always 
commensurate 0, 0, with a very unusual scaling be- 
havior of the magnetic response p,l9(. There is, however, 
a similarity between these compounds: the long-range 
Neel order is destroyed at rather close values of doping, 
x = 0.02 in LSCO and x = 0.03 in LCLioQIl^. 

It is clear that the origin of the differences between 
Sr and Li doping is in the fact that Sr substitutes La 
ions outside of the Cu02 planes while Li ions substitute 
in plane for Cu ions. Every dopant injects a hole into 
the plane, but the Sr ion seems to be a relatively weak 
hole scatterer while the Li ion interacts strongly with the 
injected holes. To understand the differences between 
these two systems we first consider the problem of an 
impurity bound state formation in a strongly correlated 
Cu02 plane. Then we generalize our solution for a finite 
concentration of holes and find that there is a critical con- 
centration x c above which a spin-shard state is formed. 
We show that the formation of this state can explain the 
observed properties of LCLiO. 

Coulomb trapping of a single hole. The problem of a 
single hole trapped by a Sr ion has been solved recently 
in Ref . [ll| • I n this work we extend that solution to the 
case of Li doping. Similarly to the case studied in Ref. 
[ll"l | we perform the analysis within the t— t'— t" — J model 
that is assumed to describe the low energy, e < At ~ leV, 
dynamics of Cu02 planes. The parameters of the effec- 



tive t—t' — t" — J model are well known: from the Raman 
scattering studies it was found that J w 125 meV; 
following the calculations of Andersen et al. |l3| we set 
t/J = 3.1, t'/J « -0.5, and t" / J » 0.3. We measure the 
energy in units of J, length in units of the lattice spacing 
a, and set % = 1 = ks- The kinetic terms t' and t" are 
small compared to t and our results for LCLiO are not 
qualitatively sensitive to their precise values [l^. Nev- 
ertheless, these matrix elements are quantitatively im- 
portant because they describe hopping within the same 
magnetic sublattice and influence the hole dispersion. 

Because of the strongly interacting environment where 
the hole moves, the hole dynamics is strongly dressed 
by quantum spin fluctuations with momenta q ~ (7T,7r). 
This dressing leads to the effective hole (quasihole or 
small magnetic polaron) that describes dynamics of the 
system at even lower energies, e < J. This part of 
the problem is well understood and is independent of 
the presence of the impurity. In momentum space, the 
minima of the quasihole dispersion appear near ko = 
(±7r/2, ±7t/2), which are the centers of the four faces of 
the reduced magnetic Brillouin zone (MBZ). In vicin- 
ity of these points the dispersion is quadratic: et ~ 
j3\k 2 /2 + 02k%/2, where k is defined with respect to ko, 
ki is perpendicular to the face of the MBZ, while k2 is 
parallel to it. In the self-consistent Born approximation 
[l5j the hole dispersion is isotropic, /3i s» fa = /3 » 2.2, 
and has a finite quasiparticle residue, Z ps 0.34. 

In the presence of a negatively charged impurity, the 
hole is trapped near the impurity ion by the Coulomb 
potential — e 2 / (t e \Jr 2 + d 2 ), where e e ~ 30 — 100 is the 
effective dielectric constant of La2Cu04 0, and d is the 
distance of the ion from the Cu02 plane. In the cases 
discussed here we can set d = since its effect is rather 
small. The most important difference between the Li and 
the Sr doping is that the hole cannot penetrate on the 
Li site because of a strong hard core repulsion. Thus, we 
can model the potentials for Sr and Li as: 



Sr : V(r) 



Li: V(r) 



+ GS(t) . 
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where G is a very large positive constant. The (5-function 
term in the Li potential enforces a node of the bound- 
state wave function at the origin. So the ground state 
wave function in the Li case can have either s-wave with 
a node, or p-wave symmetry. 

Because of the presence of the bound state the local 
magnetic configuration relaxes, leading to the formation 
of a spin texture: the spin at every site can rotate by an 
angle Qi with respect to the orientation at r — oo. There 
are two sublattices, "up" and "down", so that: 
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_ e tfl(r t )m-«r/2| 
- _t0(rj)m-ff/2| 



t), i £ "up" sublattice, 
|), j € "down" sublattice. 



Here m = (cos a, sin a, 0) with arbitrary a is the "di- 
rector" of the state. The director is orthogonal to the 
magnetization plane. Note that the directions in spin 
space are completely independent of directions in coordi- 
nate space. Because of the sublattice structure, the wave 
function of the hole tp(r) has two components, pseudospin 
up and pseudospin down, corresponding to a hole moving 
in the up and down sublattices, respectively. The total 
energy is of the form 0, 0] : 
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where p s /J ~ 0.18 is the spin stiffness of environ- 
ment and e is a unit vector orthogonal to a given face 
of the MBZ. If the hole "resides" in the pocket near 
k = (7r/2,7r/2) then e = (1/^2,1/^2), and if the 
hole "resides" in the pocket near ko = (7r/2, —tt/2) then 
e = (1/V2) — 1/V2)- There is coupling between differ- 
ent pockets, but this is an exponentially small effect and 
can be safely neglected. Finally, g = Zt « 1.05 is the 
effective hole-spin-wave coupling constant. The solution 
of is of the form 



^ (r) = 7=2 
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where, 



X(r) 
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with k and k! as variational parameters. Note that these 
solutions have s-wave symmetry. We have checked that 
the s-wave state is lower in energy than the p-wave one. 

A variation of the energy with respect to 8 leads 
to, 

V 2 9 = V2g/ Ps (e-V) X 2 (r) , 




FIG. 1: a: The radial charge density 2ivrx 2 versus radius, b: 
The angle 9 (radians) of the background spin rotation versus 
radius. The angle is given for direction r||e. The solid and the 
dashed lines correspond to Sr and Li impurities, respectively. 



which has the solutions: 

9 (e • r) 



Sr : 9 



Li : 9 



1 - e - 2Kr (l + 2«t)] 



9 (e 2 r) [l-e~ 2K ' r (l + 2^V 
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Substitution of these solutions together with (J2J and © 
in Eq. , and the minimization of energies with respect 
to k and k' leads to: 

2e 2 /(e ef 3) (3 K 2 ( g 2 \ 

27T/3 Ps ) 

, 2e 2 /(e e f3) /V 2 / 5g 2 

JjI ■ « = — : — , -C/Li = ^— I 1 — 
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The effective dielectric constant e e is not known accu- 
rately. Therefore, to determine k we rely on the esti- 
mates that follow from the hopping conductivity Q, Il5| : 
k w 0.4 w O.lA -1 . Then, from (JHJ, one finds that 
k' « 0.31 w 0.08A" 1 , £: Sr /J w -0.09 (i.e., £7 Sl « -11 
meV), and E u /J « -0.023 (i.e., £ Li w -3 meV). Al- 
though the binding energies are subject to uncertainties 
associated with the value of the dielectric constant, the 
ratio k/k' ~ 1.3 and hence relative sizes of Sr and Li 
bound states are reliable. According to measurements 
fl?| of the dielectric response of LCLiO, the charge exci- 
tation gap at x = 0.023 is A w 1 meV. This value is in 
a reasonable agreement with our estimate l-Eul ~ 3 meV 
especially having in mind that our result is relevant only 
to the very low doping, x <C 0.03. We will argue below 
that the gap decreases with doping switching from one 
regime to another at x ~ 0.03. 

Plots of radial charge density 27rrx 2 (r) for Sr and Li 
impurities are given in Fig. la. In Fig. lb we show the 
angle of the spin rotation of the antiferromagnetic back- 
ground for Sr and Li impurities (see Eqs. (J2J). Inside 
the "atomic core", r -C 1/k, the Sr solution (@J de- 
scribes a uniform (1,±1) spiral, 8 \f2gn 2 {e ■ r)/(7rp s ), 
while the Li solution Q gives a very small rotation an- 
gle, 9 ~ gK 4 r 2 (e ■ r)/p s . This result indicates that while 
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Sr impurities induce incommensurate magnetic textures, 
Li impurities do not lead to a significant change in the 
commensurate state. 

Destruction of the long-range order. According to Eq. 
(0J, at large distances, r ^> the Sr and Li ions 

generate a similar disturbance in the antiferromagnetic 
background: Sri = [no x m]8 = M(e ■ r)/(27rr 2 ), where 
n = 5 n + no is the Neel unit vector, no = n(r = oo), 
and |M| = \f2gj p s w 8 is the "dipole moment" of the 
bound state. It has been shown in Refs. that the 

dipoles frustrate the antiferromagnetism and hence lead 
to destruction of the Neel order at some critical concen- 
tration, x c . Since M is more or less the same for LSCO 
and LCLiO, we would expect x c to be the same in these 
two compounds. For M « 8 one finds x c « 0.02 [llj| . 
However, for x ~ 0.02 the average separation between 
ions is r ~ 7a and hence, according to Fig. lb the ef- 
fective dipole moment of the Li impurity is somewhat 
smaller than the one of Sr. Thus, the change in the dipole 
moment leads to a value of x c « 0.03 for LCLiO that is 
slightly larger than the value of LSCO (x c ps 0.02). This 
result explains the destruction of Neel order in LCLiO. 

In spite of similarity at very small x, the physics at 
x > x c in LSCO and LCLiO is quite different. It has been 
argued in Ref. 11] that for 0.02 < x < 0.055 LSCO re- 
mains insulating because the bound states have not per- 
colated. Nevertheless, because of the long-range dipole- 
dipole interaction between the bound states, the insu- 
lating spin-glass-like state with (1,±1) spiral ordering is 
established 20] . Due to the Dzyaloshinskii-Moriya inter- 
action the spiral direction is pinned to the orthorhombic 
b-axis. The size of the Sr-hole bound state shown in 
Fig. la is consistent with this picture. At x w 0.055 the 
Sr-hole bound states percolate and, due to the Pauli prin- 
ciple, there is a rotation of the spin spiral direction by 
45° to the (1,0) or (0,1) incommensurate order At 
x > 0.055 the charge distribution is more homogeneous. 

According to Fig. la the size of the Li bound state is 
a factor 1.5-2 larger that that for Sr. Therefore in 
LCLiO there is no window of doping between destruc- 
tion of the Neel order and percolation of bound states. 
Hence, percolation occurs at x ~ x c ~ 0.03. However, 
after the percolation the system remains extremely in- 
homogeneous because of strong scattering of holes by Li 
ions. The hole localization length should be of the order 
of the mean free path that is of the order of the distance 
between impurities. Since the localization length is the 
same as the average separation between holes, the regions 
occupied by a single hole and regions occupied by two (or 
a few) holes are equally probable. It has been shown in 
Ref. 01 that single hole regions induce (1, ±1) spirals 
while few hole regions induce (1,0) or (0,1) spirals. Thus, 
at any doping above the percolation the system remains 
completely frustrated and no collective spiral order can 
be established. The system consists of weakly and ran- 
domly coupled antiferromagnetic clusters. We call these 



clusters antiferromagnetic shards. 

This picture leads to a prediction for the nuclear 
quadrupole resonance (NQR) in LCLiO. It is known that 
in LSCO the Cu NQR A-line is shifted with respect to 
its position in the parent compound La2Cu04, and the 
shift is proportional to doping j^]. The A-line is wide, 
reflecting a non- uniform charge distribution In the 
spin-shard state of LCLiO (i.e. at x > 0.03) one should 
expect two distinct NQR lines of comparable intensity 
(two Cu lines and two O-lines). The first one originates 
from antiferromagnetic shards. It is likely that the fre- 
quency of this line is close to that in the parent com- 
pound. The second line originates from regions occupied 
by holes and, therefore it should have a higher frequency 
and larger width. We note that inhomogeneities have 
already been observed in La-NQR in the ordered phase 
(x < 0.025) [2^]. We assign those to the formation of 
bound states discussed here. 

The spin-shard state. Below percolation, x < x c the 
system retains the long-range antiferromagnetic order. 
Above percolation, x > x c there are only finite antifer- 
romagnetic shards with N spins because all exchange in- 
teractions are short range. In fact, there is a distribution 
of shard sizes and the probability of finding a shard of 
size N is given by percolation theory jiij : 

P(N, x) oc exp{~N/N(x)} , (6) 

where N(x) is the average shard size at doping 
x. Because of magnetic anisotropics in the system 
(Dzyaloshinskii-Moriya and crystal field) the low energy 
physics of the shard can be thought as a two-level sys- 
tem with a tunneling splitting, A (AT), given by the WKB 
expression: 

A(N) oc ujq exp{-77V} , (7) 

where ujq oc J is an attempting frequency and 7 is a 
parameter that depends on the magnetic anisotropics. 
The combination of © and J7J) shows that there is a 
distribution of tunneling splittings in the problem that 
has the form: 

P(A) oc A"" 1 , 

where a = 1/(7 AT (x)) is a non-universal quantity which 
depends on doping. The neutron scattering is peaked 
at q = (tt, it), and the imaginary part of the average 
susceptibility is given by: 

X >, B, T) = J dAP(A) x ^ shard (c, B, T, A) , (8) 

where fi = 0, 1 refers respectively to longitudinal and 
transverse polarizations of the neutron magnetic field 
with respect to the applied uniform external magnetic 
field B. The shard susceptibility is given by the two-level 
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system expression (u> > 0): 
(A 2 ) 1 ~^(7W 2 )^ 

X^shard K A 2 + J^2 

tanh(^^)^-VA^T^), 

where .M(-B) oc ^lbVNB is the Zeeman energy (since a 
shard of ./V spins is a random system, there are on average 
^/N more spins in one sublattice than in the other leading 
to a net shard moment). Hence, the integration in 
gives the following result: 

X>,B,T) = ^—F^ U /T,MM , (9) 

where 

F^x, y) ex x Q -y^(l - y 2 )"/ 2 ^ tanh(ar/2)e(l - y) 

is a scaling function [23, where Q(x) = 1(0) if x > 
(x < 0) is the Heaviside step function. In a unpolar- 
ized neutron scattering experiment both longitudinal and 
transverse responses appear. When y ~ 1, a s in the case 
of nuclear magnetic resonance (NMR) [26j, the trans- 
verse response is the dominant one. In this case, since 
we expect a < 1, the transverse susceptibility (/1 = 1) be- 
comes singular at y = 1 while the longitudinal response 
(/j = 0) vanishes. This resonant effect should be easily 
observable in neutron scattering experiments in a finite 
magnetic field because the damping to the shard motion 
should be weak due to the insulating nature of the back- 
ground. Equation @ is in agreement with zero magnetic 
field (fi = 0, y — > 0) experimental data of Ref. |8| with 
a w 0.35 for x = 0.04 and a w 0.06 for x = 0.06 and 
x = 0.1. By performing the same experiments in a mag- 
netic field one can test our prediction for (JSJ. 

In reality the antiferromagnetic shards are not decou- 
pled. Shards from different CuC>2 planes interact because 
of the interplane exchange J± . Due to this interaction the 
system must freeze to a disordered 3D spin-shard glass 
at the characteristic temperature T30 ~ J±N ov , where 
Nov ~ N(x) is the number of spins in two shards that 
overlap along the c-axis. Thus, Eq. JJjJ is relevant to 
the spin-shard state at T > T 3 d where shards fluctuate 
independently. At T < T 3 £> the shards freeze. However, 
Eq. © is still valid, but instead of the current tem- 
perature T one must substitute the freezing temperature 
T 3D - Thus, for T < T 3D one finds that x"(u,T) should 
be roughly temperature independent, in agreement with 
the experimental data 0,0. 

In conclusion, based on an analysis of the extended 
t — J model, we propose the theory of antiferromagnetic- 
spin-shard state for La2Cui_ a; Li, c 04. The theory ex- 
plains why the long-range Neel order is destroyed at 
x c = 0.03. At x < x c physics of the system is driven 
by Li-hole bound states. At x = x c the Neel order is de- 
stroyed and simultaneously the bound states percolate. 



At x > x c the holes in L^Cui-xLi^C^ remain localized 
breaking the global antiferromagnetic state into array of 
weakly coupled 2D antiferromagnetic shards. Due to the 
interlayer interaction the shard state freezes to a glassy 
state (spin-shard-glass) below some characteristic tem- 
perature T30. While our theory describes quite well the 
available experimental data, it also predicts two indepen- 
dent effects: the existence of two distinct Cu(0)-NQR 
lines, and a magnetic field scaling of the neutron scatter- 
ing data. 
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